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Abstract 

The strong dependence of the mass variance parameter, ag, on the adopted cluster 
mass-temperature relation is explored. A recently compiled X-ray cluster catalog 
and various mass-temperature relations are used to derive the corresponding values 
of as- Calculations of the power spectrum of the CMB anisotropy induced by the 
Sunyaev-Zeldovich effect and cluster number counts are carried out in order to assess 
the need for a consistent choice of the mass-temperature scaling and the parameter 
(Tg. We find that the consequences of inconsistent choice of the mass-temperature 
relation and as could be quite substantial, including a considerable mis-estimation 
of the magnitude of the power spectrum and cluster number counts. Our results can 
partly explain the large scatter between published estimates of the power spectrum 
and number counts. We also show that the range of values of the power-law index in 
the scaling Ci ~ Cg"'' deduced in previous studies is likely overestimated; we obtain 
a more moderate dependence, ~ (t|. 
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1 Introduction 



Calculations of the power spectrum of CMB anisotropy induced by the Sunyaev- 
Zeldovich (S-Z) effect, and of cluster number counts, involve various intrinsic 
cluster quantities, as well as cosmological and cluster mass function parame- 
ters. Specifically, the level of S-Z signal in a cluster is directly related to the 
intraclustcr (IC) gas temperature and density distributions. The cosmological 
and large scale structure parameters enter in the expression for the collec- 
tive mass function of clusters. Most commonly, this function is assumed to 
have the Press & Schechter (1974) form, but other, more physically estab- 
lished mass functions, are also widely adopted. The mass function is usually 
normalized by specifying an observationally deduced value of Ug, the rms den- 
sity fluctuation on a scale of Sh~^ Mpc. The value of erg may be deduced in 
several different ways, perhaps most directly from a comparison between the 
observed two-point correlation function of galaxies or clusters of galaxies, and 
the spectrum of the primordial density fluctuation field on the relevant scales, 
recently employed in the analysis of the large SDSS sample (e.g. , Tegmark et 
al. 2003). Measurements of CMB anisotropy yield a globally averaged value 
of (Tg, most recently deduced from the first year WMAP measurements (Ben- 
nett et al., 2003). Additionally, erg has been deduced from fitting X-ray cluster 
data (luminosity or temperature functions) to the corresponding theoretical 
models, with the latter derived from a mass function upon assuming a specific 
mass-luminosity or mass-temperature relation. 

Deduced values of cg clearly depend on data quality, modeling, and analy- 
sis methods, as is reflected in the current wide range, 0.5 ^ o"g < 1.2. Joint 
analysis of the SDSS and WMAP databases yields the current 'best' (with 6 
free 'vanilla' parameters) value, cxg = 0.917iQ;o72 (Tegmark et al. 2003). This 
is roughly in agreement (e.g., Wu 2001) or somewhat lower (Viana & Liddle 
1999) than results from X-ray cluster studies. The wide range of values of 
Ug has signiflcant implications on the level of various cosmological quantities. 
Here we explore some of its consequences for the abundance of clusters, and in 
particular the power of S-Z anisotropy and cluster number counts, quantities 
that depend steeply on ug. Since we concentrate here on clusters, the emphasis 
will be on values of erg deduced from cluster X-ray measurements, particularly 
of gas temperatures. Pertinent key features of such analyses are properties 
of the cluster sample (e.g. , completeness, temperature and rcdshift ranges), 
the adopted mass function model, and the mass-temperature relation used to 
convert the mass function into a temperature function. By comparing the lat- 
ter function to the observed temperature function, (jg (and other parameters, 
such as D,m) can be determined. 

In this paper we focus on consequences of the choice of the temperature-mass 
(TM) relation on (jg, and their impact on the predicted S-Z angular power 
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spectrum and cluster number counts. A quantitative assessment is given of 
the impact of inconsistent choice of the TM relation and (Tg on the predicted 
levels of the the S-Z anisotropy and the number of clusters to be observed by 
the Planck satellite. Our main goal here is to improve some of the technical 
aspects in the calculations of these quantities, and to partly explain the large 
variance among previous estimates of the S-Z power spectrum. 

The paper is arranged as follows: the need for a consistent choice of as and the 
TM relation, the method employed to extract as from the observed tempera- 
ture function, and the variation of as with the choice of the TM relation are 

described in § 2. The importance of consistent choice is addressed once more 
in § 3 and § 4, where its consequences on the S-Z angular power spectrum and 
cluster number counts are quantified and further elaborated upon in § 5. 



2 TM-as Relation 



2.1 Method 



The calculation of the S-Z power spectrum involves the basic cosmological 
and large structure parameters, properties of IC gas, and the cluster mass 
function (MF). Of the various methods that have been used to calculate the S- 
Z anisotropy, we employ the approach described in Colafrancesco et al. (1994, 
1997), which we have recently updated and generahzed (Sadeh & Rephaeli 
2003). The calculation will not be described here; full details are given in the 
latter paper. 

In a cluster temperature-mass relation (TMR) the gas temperature is a mono- 
tonically increasing function of mass; T oc M^/^ in virial equilibrium. The 
temperature function may be written in terms of the MF by 

dT 

N{T) dT = N[M{T)] — dM. (1) 



Since at the high mass end the MF is to a good approximation a monotonically 
increasing function of as, it can be formally shown (Sadeh & Rephaeli 2003) 
that: 

1. TM relations which predict higher (lower) temperatures for a given mass 
yield lower (higher) as, and therefore, lower (higher) cluster density. 

2. Calculations of the S-Z power spectrum employing a specific TMR, and a 
value of as different from the one that would have been extracted by adopting 
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this relation, will result in mis-cstimation of S-Z power. The power level will 
be overestimated (underestimated) when, for a given TMR, the adopted value 
of (Tg is higher (lower) than the value corresponding to that particular TMR. 

We adopt the procedure developed by Pierpaoh et al. (2001) for the deter- 
mination of as from an observationally based temperature function. Given a 
list of X-ray cluster temperatures and their corresponding redshifts, the tem- 
perature function can be deduced. The theoretical MF - with erg as the free 
parameter - is converted into a temperature function given an assumed TM 
relation. This is calculated at redshift z — 0.057, the median redshift of the 
particular X-ray catalog used. The temperature range in the catalog is divided 
into a large number of bins, and the mean occupation number of clusters in 
each of these bins is computed by multiplying n(T) dT at the bin temperature 
by the effective volume occupied at this temperature, i.e., the volume of space 
to which clusters in the bin can be observed, denoted by . A large number 
of mock catalogs is now devised by using the listed temperatures as the mean 
of a normally distributed variable, and their la errors as its variance. The 
temperature bins are now scanned so as to determine whether any of these is 
occupied by clusters in the mock catalog. Occupation numbers are assigned to 
each one of the bins. Assuming a Poissonian distribution of the cluster occupa- 
tion numbers, the probability of finding Ni clusters in bin i is estimated using 
the Poissonian distribution, with n{T) dT Vi interpreted as the distribution 
mean. Finally, the likelihood function of obtaining the 'observed' occupations 
in all bins simultaneously is constructed, and maximized against as- 



2.2 TM Relations 



We employ five TM relations found in the literature; three of these have the 
same functional form, differing only in the temperature normalization for a 
fiducial cluster at redshift z — with mass 10^^ h"^ M©. The five functions 
are identified here as models 1-5: 

Model 1 (Colaf ranees CO et al 1997): 

( M 



A(l]o,z) 
A(Qo = 1,-^ = 0) 



I/O 



(2) 



Model 2 (Tomita 2003): 



r = 5.27(l + ^) 



M 



2/3 



1015 /i-i Mq 



1/3 



Pvir (^) 

187r2 



1/3 



(3) 
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where 7 = 1.2. 

Model 3 (Henry 2000): 



_ 7.98 A;e1/, 
r=— ^ (1 + ^) 



M 



2/3 



TM 



1015 h-^MQ^ 



1/3 



187r2 



nl/3 



(4) 



where Ptm — 1-21. 



Model 4 (Pierpaoli et al. 2001): 



M 



2/3 



10^5 /l-l Mr 



0/ 



(A,£;2)i/3 _ 2 



(5) 



where E'^ = Q^(l + 2;)^ + Qa- 

Model 5 (Molnar & Birkinshaw 2000): 



r = 2.76/3 



2/3 



27r 



2/3 



sinhr] — r] 



+ 



Up + 3 



, (6) 



where (3 and Up are taken to be 1 and —1.4, respectively, (n^ is the 'effec- 

r 2 ~ 

tive' power index on cluster scales), rj — cosh~^ 



Qoil + z) 



ilm{z) - 1 



, and ^m{z) — 



no{i + z) + {i-no)' 

The first three relations may all be written in the form 



T^a{l + z) 



M 



2/3 



1015 h-^MQ^ 



1/3 



187r2 



1/3 



(7) 



where a — 5.78,6.24,6.59 in models (1),(2), (3), respectively. In a fiat uni- 
verse with A = 0, a is equal to the temperature of a local cluster with mass 
10^5 /i^i Mq. In the currently fashionable ('standard') ACDM model, substi- 
tuting z = 0.053 in the above relations, the temperatures associated with 
clusters having this mass are 3.43 fceV, 3.71 /ceV^, 3.92 fceV, 6.19 fceV^, 8.49 fceV^, 
respectively. 
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2.3 Results for as 



To verify the validity of the first assertion made in the previous subsection, 
(jg has been estimated using equation 7 with 5.5 < a < 7.0. For each value 
of a ten runs were made and average values of (Tg were computed. A more 
statistically complete analysis would require a much larger number of runs for 
each value of a, but the consistency of the results suggests that this is not 
needed for our purposes here. The input and output of the runs are tabulated 
in table 1 and illustrated in figure 1. It is clear from the plot that (as expected) 
(Tg is quite sensitive to the choice of the TM relation. In fact, a relative vari- 
ance of ~ 30% in the temperature scaling - which is in accord with typical 
observational uncertainties - induces a relative change of ~ 20% in the value 
of (Tg. This has important practical consequences given that the mass function 
is extremely sensitive to this parameter. Moreover, these results are in accor- 
dance with the qualitative argument summarized above: A TM relation which 
associates a higher temperature with a given mass yields a lower value of cig, 
and vice versa. The calculated five values of cxg are 1.32, 1.23, 1.18, 1.00, 0.68, 
for models 1-5, respectively. 



3 Dependence of the Power Spectrum on erg and the TMR 

The S-Z power spectrum was calculated (Sadeh & Rephaeh 2003) employing 
the five models described above. We first calculated the power for each of 
the five TM relations and its matching ag, and then crossed these inputs by 
taking the TM relation of model i and o"g corresponding to model j, where 
i ^ j = 1 — 5. Results of the first calculations are shown in figure 2; these 
reflect the tight relation between the magnitude of the power spectrum and 
the cluster population. Note that although lower temperatures (implied by a 
'colder' TM calibration) will tend to reduce the power, the influence of the 
higher erg dominates, resulting in a higher level of power. 

The results of the second set of 'cross' calculations are described in flgure 3. 
Each panel contains three curves; the solid curves represent calculations em- 
ploying a specific TMR and its corresponding cTg. The dash-dotted curves show 
results of the calculations employing the same TMR as for the solid curves, 
but a different cxg. Finally, the dashed curves show results for the same cxg as in 
the dash-dotted curves, but with the TMR that corresponds to this parameter, 
thereby providing a 'correction' to the results illustrated by the dash-dotted 
curves. In order to further clarify the results, the first row will be considered 
as an example: In the first panel (2nd from the left) the TMR of model 1 is 
used; here the temperatures are lowest, and erg assumes the highest value; the 
solid curve describes this case. Now erg is changed to the value extracted using 
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Table 1. Values of < crs > for the adopted TM relation 



model 2, i.e., a lower value. However, in combination with the TMR of model 
1, the power spectrum is underestimated, since the TMR of model 2 should 
have been used, rather than that of model 1. Indeed, model 2 predicts higher 
temperatures than model 1, and therefore the correct power spectrum is rep- 
resented by the dashed curve. Obviously, as the difference between the correct 
(78 and the one employed in the calculation increases, the degree of underes- 
timation (or overestimation) also increases, as can be seen in the next panels 
of the first row. The relative errors introduced by the inconsistent choice of 
these two elements can be readily calculated; this is presented in figure 4. As 
anticipated, the panels lying on both sides of the diagonal show lower relative 
errors than those located further away from the diagonal, where differences 
between adopted and correct values of are larger. Errors of up to a factor 
of ~ 5 are caused by this parameter mismatch. The relative errors among the 
first three models are constant along the multipole axis due to their functional 
uniformity. In other cases the relative error changes with £. 

The dependence of the S-Z power on cig has already been discussed in previous 
works, it is usually demonstrated by arbitrarily changing its value. Although 
the change is not great (as the ones required in order to induce a relative 
error of hundreds of percents as seen in figure 4), a difference of < 20% (e.g. 
(78 = 0.9 — 1.1) or more is common. Examples for large errors induced by a 
change of ~ 18% in (jg may be seen in the third panels of the third and fourth 
rows of figure 4: in the third row the correct and unmatched values of are 
1.18 and 1.0, respectively, and a relative error of ~ 40 — 50% appears around 
I — 1000. In the fourth row the correct and unmatched values are 1.0 and 
1.18, respectively, and a relative error of > 50% is induced. 

The steep dependence of the S-Z power on ag was quantified by the proportion- 
ality Ci oc cTg~^ (e.g., Komatsu & Tetsu 1999). Our results indicate that this 
dependence is somewhat more moderate, which is to be expected since an in- 
crease of (78 implies a decrease in temperature by virtue of the anti-correlation 
between these two parameters, as has been shown above and illustrated in fig- 
ure (1). For the three representative value of the multipole (marking the range 
in which the peak of the S-Z power spectrum occurs) £ = 1000, 2000, 3000, 
we plotted the corresponding C(S in models 1-5 against their corresponding 
(Tg. Best fits describing the three sets of five data points all have slopes ~ 4, 
smaller than the values (~ 6 — 7) previously reported. This is illustrated in 
figure (5). For comparison, we include figure (6), in which the same analysis 
is applied, but here the TMR models 2-5 are used in conjunction with the cxgs 
extracted using models 1-4. The corresponding power spectra are represented 
by the dash-dotted curves in the four panels under the central diagonal of fig- 
ure (3). Recall that the resulting power spectrum in this case is overestimated 
due to the fact that the TMR used in its evaluation would necessitate a lower 
(Tg than the one used in the calculation. Thus, a steeper slope is expected 
in the logQ — log (78 plot, which is indeed the case. The resulting slopes are 
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actually steeper (~ 13) than a ~ 6 — 7, but the essential result of this analysis 
is that the scaling — as is indeed sensitive to the adopted TMR and its 
normalization. 



4 Dependence of Cluster Counts on ag and the TMR 



As is well known, the fact that the S-Z effect is independent of redshift makes 
it an (almost) ideal tool for the detection of distant clusters. Cluster surveys 
are planned by several S-Z projects, culminating (according to current plans) 
with the full sky mapping by the PLANCK sateUite. The 143 GHz and 353 
GHz channels of the HFI experiment on PLANCK are particularly suitable 
for detecting large numbers of clusters. The S-Z flux from a cluster is given 
by (e.g. Colafrancesco et al. 1997) 

AF,(7,) = I R,{\^-ielaB)AI,{^)dn, (8) 



where A/j, is the change of the CMB spectral intensity due to the thermal 
component of the S-Z effect. In the non-relativistic limit, AI^, = g{x) = 
x'^e^{e^ — 1)^^ coth(a;/2) — 4], where Tq is the CMB temperature, and y 
is the Comptonization parameter integrated along a los. The more accurate 
expression for AIi, includes relativistic corrections (Rephaeh 1995, Itoh et al. 
1998, Shimon & Rephaeli 2003). Since our main aim here is to display the de- 
pendence of cluster number counts on erg and the TMR - rather than precisely 
predicting the expected number of clusters - we use the simpler non-relativistic 
limit of A/j,. We can then write for the S-Z flux 

AF, = ^^^^g{x)yo J Rs{\^ - iel as) ■ Cillil z) dil, (9) 



where 



0T(w 



tan" 



P 



1 - jO/pOc) 

\ 1 + {0/0cy 



(10) 



The measured signal is actually weighted over the spectral response of the 
detector, Ei,, usually approximated by a Gaussian, G^- Therefore, 
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For simplicity, Gy is taken here to be uniform over the passband centered on 
the central frequency. For the PLANCK/HFl 143 GHz channel, beam size is 
8' and lS.v jv = 0.25. The predicted number of clusters with flux greater than 
AFi, can now be calculated using 

N{AF^) = jr^^dz j B{M, z)N{M, z)dM. (12) 



The lower limit of the mass function integral is such as to correspond to the 
limiting flux from a cluster with mass M, which we take to be 30mJy, the 
flux limit reachable in a one year observation. To evaluate this limiting mass, a 
binary function B{M, z) is introduced. B{M, z) = 1 if the flux corresponding 
to a given mass (and redshift) is greater than the limiting flux; otherwise, 
B{M,z) =0. 

The predicted number of detectable clusters with mass M at redshift z is 
proportional to the MF evaluated at these parameters. This depends on the 
integrated Comptonization parameter along all lines of sight to the cluster. 
Consequently, the arguments presented in § 2 are relevant also in this case. In 
figure 7 the predicted cumulative number counts in models 1-5 are shown. A 
strong relation between the MF normalization and the cluster counts is clearly 
manifested. It may also be inferred from the figure that the largest contribution 
to the counts comes from clusters lying at redshifts < 0.4. Since the redshift 
distribution of clusters cannot be determined from S-Z measurements, the 
most relevant quantity is indeed the cumulative cluster counts at low redshifts. 
Results of the 'cross' calculations are presented in figure 8. The curves in 
figure 8 resemble those in figure 3 in the sense that the differences between the 
values represented by the 'mismatched' and 'correct' curves increase farther 
away from the diagonal. The relative differences between the properly and 
improperly calculated numbers span a range of ~ 10 — 70% at the lowest 
redshifts, corresponding to the maximal cluster counts that would be detected 
in an all-sky S-Z survey. 

Although not observable, the larger differences exhibited at higher redshift 
refiect the higher sensitivity of the MF to Ug at the high-mass end. At such 
high redshifts only the most massive clusters would be detected, owing to 
their larger fiuxes, and therefore any change in the normalization that does 
not match the TMR used to evaluate the fiux, would result in larger devia- 
tions of the counts with respect to a model employing a specific TMR and 
its corresponding normalization. A comparison between the results illustrated 
in figures 3 and 8 suggests that an inconsistent choice of as and the TMR 
results in a rather pronounced over- or under-estimation of the power spec- 
trum, and a milder (but significant) mis-estimation in the cumulative number 
counts, in particular at the lowest redshifts. This difference is clearly due to 
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the steeper (gas) temperature dependence of the power (oc T^) as compared 
to the corresponding hnear dependence of the flux. 



5 Conclusion 

Predicted profiles and levels of the angular power spectrum of the S-Z effect 
and of cluster number counts are strongly affected by modeling details and the 
particular choice of parameter values which explain the large variance in the 
values of these quantities in the literature. In this paper we focused on the tight 
dependence between the TMR and erg. We demonstrated that selecting these 
independently is generally not self-consistent and leads to considerable mis- 
estimation of the power spectrum and cluster counts. This can be is avoided 
when as is determined from CMB and large scale structure surveys that do 
not involve cluster TMRs. 
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Foundation grant 729\00-03 at Tel- Aviv University. 
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Fig. 1. (Js as function of the parameter a. Each data point represents an average 
over ten runs. The fitted curve is plotted only for visual purposes. 
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Fig. 2. The angular power spectrum of the S-Z effect calculated using the five 
M-T relations of models 1-5, and their matching ags as presented in § 3. Solid, 
long dashed, short dashed, dashed-dotted and dotted curves correspond to models 
1,2,3,4 and 5, respectively. 
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Fig. 3. The S-Z power spectrum calculated using unmatched TM relations and as 
normalization. Units of the abscissa and ordinate are £ and 10~^^ £{£+ 1) C^/27r, re- 
spectively. Rows 1-5 correspond to models 1,2,3,4,5. Columns 1-5 correspond to the 
correct ag for the model. Solid curves depict the 'diagonal' terms, that is, calcula- 
tions carried out with a specific TMR and its corresponding erg. Dash-dotted curves 
describe calculations made with the same TMR, but unmatched as- Finally, the 
dashed curves illustrate the 'corrected' power spectrum, i.e. keeping the unmatched 
as but now using the matching TMR. 
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Fig. 4. As in figure (3), but here the relative differences between the dashed and 
dash-dotted curves (i.e., consistent vs. inconsistent calculations) are specified in 
percentage units. 
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Fig. 5. The dependence of Q on as, illustrated in a Log-Log plot. Asterices, x-marks 
and plus symbols represent Qs at five different ags, corresponding to models 1-5, 

for i = 1000, 2000, 3000. The dashed, solid, and dash-dotted curves depict the best 
linear fits for each of the 5 data point sets. The two short lines appearing at the 
bottom of the plot are of slopes 4 and 6, as indicated in the figure. 
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Fig. 6. The same as figure (5), but here the TMR models 2-5 are used in conjunc- 
tion with non-matching ags, extracted from models 1-4, respectively. Thus, both 
higher values of a% and temperatures scaled from the TMR are used, leading to an 
overestimation of the power spectrum. 
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Fig. 7. S-Z cluster number counts calculated using the five TMR described above, 

and their matching values of as- As in figure (2), solid, long dashed, short dashed, 
dash-dotted and dotted curves correspond to models 1,2,3,4 and 5, respectively. 
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